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Abstract 

A method for finding an optimum «-dimensional commutative group code of a given order M is 
presented. The approach explores the structure of lattices related to these codes and provides a 
significant reduction in the number of non-isometric cases to be analyzed. The classical factor- 
ization of matrices into Hermite and Smith normal forms and also basis reduction of lattices are 
used to characterize isometric commutative group codes. Several examples of optimum commu- 
tative group codes are also presented. 
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1. Introduction 

The design of spherical codes for signal transmission through a Gaussian channel is a clas- 
sical problem in coding theory, where group codes have proved useful [1 2 1 since their ap- 
pearance in the pioneering work of Slepian Q. The special attention devoted to these codes is 
largely due to their symmetry and homogeneity which arise from their special algebraic structure 
|4][5]. The interest in such group codes has persisted with various studies have been developed 
10 17] [8] [9] [TO] [TT], including some proposing applications in turbo concatenated and low density 
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schemes IfTSl [131 [141 [T5l [T^l . Recently it has been shown that the Shannon capacity of cer- 
tain important channels, as the AWGN channel with m-PSK modulation, can be achieved using 
commutative group codes IfTTl and they will be focused here. 

One of the underlying difficulties in the design of a group code is the finding of an ini- 
tial vector which maximizes the minimum distance of the associated code, for a fixed group of 
orthogonal matrices; the so called initial vector problem. This problem still does not have a 
general solution, although various important cases have been studied, including reflexion group 
codes |[T8l and permutation group codes [19|. Besides, Biglieri and Elia have shown ll20ll that 
for cyclic group codes the problem can be formulated as a linear programming problem. Here 
we extend their ideas and show that for any commutative group code, the initial vector problem 
can also be solved in the same way. 

Furthermore, this paper deals with the more general problem of determining an optimum 
commutative group code in R" for a given order M. We derive a two-step algorithm which 
leads to the finding of a code with maximum minimum distance for a fixed number of points M. 
Our approach explores the connection between even dimensional commutative group codes and 
lattices related to them in the half of the dimension ||2T1|2"2"1 . Using basis reduction of lattices and 
the classical factorizations of matrices into Hermite and Smith normal forms, we characterize 
a set of relevant cases to be analyzed, after discarding isometric codes. The reduction process 
presented here can also be used in the solution of other problems where lattices fl23l , in particular 
orthogonal sub-lattices, are involved; including coding and decoding process l24l l25l l26l l27l . 
image compression l28l . spherical codes on torus layers E9ll and also the enticing lattice based 
cryptography lf30l l3ll l32l . 

This paper is organized as follows. Commutative group codes and some of their properties 
are presented in Section 2. We then discuss the initial vector problem for those codes and charac- 
terize it as a linear programming problem in Section 3. The main results are presented in Section 



4, where we prove a simple, but useful, extended Hermite normal form (theorem 4. 1 1 which 
allows the characterization of isometric lattices by coordinate permutation; in this section we 



also derive theorems 4.2 and 4.5 which provide a significant reduction in the number of codes 
to be checked in the search for an optimum one. Our method is presented as a pseudo-code 
(Algorithm 1), and some examples of optimal codes in several dimensions are given. 
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2. Commutative group codes 

Let 0„ be the multiplicative group of orthogonal matrices n x n and Q n (M) be the set of all 
order M commutative subgroups in 0„. 

A commutative group code C is a set of M vectors which is the orbit of an initial vector xq on 
the unit sphere S n_1 c R" by a given G e Q n (M), i.e. 

C := Gx = {gx ,g e G) . 

We assume that C is substantial, i.e., not contained in a hyperplane. 
The minimum distance in C is defined as: 

d := min ||jc — y||= min ||gjjr — jc||, 

.v.yeC g;£C 
x i= y ft * I„ 

where ||.|| and /„ denote the standard Euclidean norm and the identity matrix of order n respec- 
tively. 

In what follows, C(M, n, d) denotes a code C in M" with M points and minimum distance 
equal to d. A C(M, n, d) is said to be optimum if d is the largest minimum distance for a fixed M 
and n. 

As is well known, the minimum distance of a group code C, generated by a finite group G, 
may vary significantly depending on the choice of the initial vector xq. Therefore, the search for 
an optimum n-dimensional commutative group code with M points requires the consideration of 
all G e Q n {M) and solution of the initial vector problem for each G. 

A well known real-irreducible representation of a finite commutative group of orthogonal 
matrices G can be stated as follows: 

Theorem 2.1. ( H33V Theorem 12.1) Every commutative group G e @ n (M) can be carried by the 
same real orthogonal transformation q into a pseudo-diagonal form: 



qgil' = ■ ■ ■ ,Rk(i),l*(t)2k+n ■ ■ ■ ,fj-(i)n\nxn, 



bij € Z, < b u ^ M and = ±1, 1 = 2k + 1, ...,n, j= 1, ... ,k, Vg; e G. 

3. The initial vector problem 

In this section we consider, for each group G e Q n (M), the search for a vector x in S"~ l 
which maximizes the minimum distance between two points in C — Gx, i.e., the search for an x 
that solves: 



max min \\giX - x\\ 

xeS»-< \gieG,gi^I„ 

This initial vector problem has been solved only in certain special cases. Biglieri and Elia have 
shown in 1 20 1 that, for cyclic groups, this search can be reduced to a linear programming problem 
(LP). Here, we extend their ideas and present an alternative formulation which also allows the 
reduction of the initial vector problem to a LP for any commutative group code. 



According to Theorem 2.1 we have: 



\\ gi x-x\\ 2 =2-2 



Considering 



we obtain 



k n-k 

Z I 1 ~ 2 sin2( M bij) ) + 4) + Z mm*) 
\j=i * j=k+\ 



y.i 



x lj-i +x li ,Hj = h-..,k 



l 2j-l T A 2j 

x 2 

X j+k 



if j = k + I, . . . ,n — k 



\\giX-x\\ =2-2 



' k n-k \ 

-2sin 2 (^^ 7 )jy /+ J] fi(i)jyj 
\j=i ' " j=k+\ 



Thus, max min||p,x - x\\ is equivalent to 

jceS"- 1 \gi*I„ 



max min 



( k 



2-2 



n-k \\ 

Z^ 281 " 2 ^^* Z Mm 

Kj=\ ' j=k+\ 



n-k 



subject to /m = l,y ; ' > 0. 



This max min problem, which is linear in y, can be reduced to the following linear programming 
problem: 



max z, 



subject to 



/ k n-k 



z <2 - 2 £ (l - 2sin 2 (£i y )]y,- + J] p{i) jyj 




n-k 



k=\ 



yi>0 



Therefore, the initial vector problem for commutative group codes is equivalent to a linear 
programming problem with n-k+1 variables. Due to the symmetry of the function sin 2 (x), in 
the case where the group G is free of 2 x 2 reflection blocks (n = 2k), the number of constraints 
can be reduced to + l)- 



4. Optimum commutative group codes 

In this section we consider a more general problem of finding a commutative group code of 
order M in W which has the largest minimum distance. To do this, we must consider all com- 
mutative groups G e Q n (M) with the respective best initial vectors and compare the minimum 
distances of the correspondent codes. 

Let us start by estimating the number of commutative group codes to be checked in order to 
find an optimum one. 

As usual, we say that two groups G and H are equivalent if they are conjugate, i.e., 



Although the set @„(M) is infinite, conjugate groups generate isometric codes. Specifically, 
given an initial vector x, G e Q n (M) and p e 0„, the group code generated by G is isometric to 
the group code generated by H = pG p', with initial vector p x. In fact, for each hi = pgip' € H, 
we must have 




G ~ H <==> 3 p eO„ ;H = pGp'. 
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i(px)-px\\ = \\(pgip')(px) - px\\ = \\pgjX-px\\ = \\gjX-x\\ 



Thus, the search for optimal commutative group codes can be restricted to groups which are 
distinct up to conjugacy. In other words, it is sufficient to consider just one representative for 



each class of the quotient Q„(M)j «, resulting in a finite set. In fact, by Theorem 2.1 for each 
G G @n{M) there exists H = qGq' in a pseudo-diagonal form i.e., for each class in quotient 
Q n (M)l «, there is a representative in the pseudo-diagonal form. Therefore, in the search for 
optimum commutative group codes, it is sufficient to consider only the set of commutative groups 
such that their matrices are in the form ([TJ. Let us denote this set by B n . The cardinality of B„ is 
clearly finite, since bjj ^ M. 

However, the set B„ still has equivalent groups and can be reduced. For instance, let G e B n 
be a group of matrices free of 2 x 2 reflection blocks, i.e., the elements in G have only 2x2 
rotation matrices as diagonal blocks. 

Let 



Gtj = 



be the j-th block of the j'-th generator of G. The block G, ; is a rotation by an ang leof febijlM). 
Note that the rotation block corresponding to M - bjj is a conjugate of the block associated with 



CO H M ) '""I M ) 

. /2x(M-b li )\ /2,r(M-6 ; )\ 
Sm l M ) COS \ M ) 



cos (t) sin (T) 
- sin (T) cos (t) 



[ -1 J G U [ -1 J • 



Therefore, up to conjugacy, we can consider b,j < M/2 in ([T]). 

Moreover, the permutation of two consecutive blocks G,y and G(,- + iy (and hence any two 

rotation blocks) results also in a conjugacy in 02k- We next consider only the set B n and also 

discard equivalent groups, as described above. 

In 11201 . Biglieri and Elia present the estimation ( ^ | for the number of cyclic groups which 

\n/2) 

must be checked in order to find an optimum one. By discarding isometric codes, as stated above, 
and also considering the Adams ' condition [34], presented next, we can give a lower estimate for 
the number of cases to be tested in the search for an optimum cyclic group code. 



6 



Adams' condition: for a fixed M and a, b e Z A we say that a and b are 
Adam's- equivalent and denote by a - b iff there exists a invertible in %m 
such that a = ab mod M. 

A generator matrix of a cyclic group G e Ozk can be defined by a vector b = {b\, b%, ■ ■ ■ ,bk) 
with < bi < M and gcd(£>i,£>2> ■ ■ ■ ,bk,M) — 1 to represent the rotation blocks. The Adam's 
relation, a - b, implies that two pseudo-diagonal matrices ([T]i with parameters defined by a and 
b generate the same cyclic group. 

(M/2\ 

Thus, the number of distinct cyclic groups is clearly less than and depends on the 

\n/2) 

number of invertible elements in Zjf, which is given by the Euler phi function of M, ip(M), 
Moreover, as pointed out above, we can restrict our search to vectors b = [b\,bz, •■ ■ ,bk), < 
bi < M/2. Based on these arguments, we can estimate the number of cyclic group codes, up 
to symmetry, by (M/2) k /w(M), which is lower than the number \ |, given in [20]. Table 1 

W 2 / i_ 

shows a comparison of these values for k = 2 and several values of M. The final column refers 
to the number of cyclic group codes effectively tested by Algorithm 1 (derived in Section 4), 
which discards additional isometric groups in order to find an optimum code. 

Table 1 : Different estimations for the number of distinct (non-isometric) order M cyclic group codes in R 4 and cases 
effectively tested by Algorithm 1. 



M 


(M/2\ 


M 2 


Algorithm 1 


\n/2) 


4ip(M) 


32 


120 


16 


14 


64 


496 


32 


26 


128 


2016 


64 


50 


256 


8128 


128 


98 


512 


32640 


256 


194 


1024 


130816 


512 


386 



In what follows, we will focus our attention on the class of commutative group codes, with 
generator matrices are free of 2 x 2 reflection blocks. Moreover, it is sufficient to consider 
commutative group codes in even dimensions, because, as pointed out in [21 1, a commutative 
group code in odd dimension, n = 2&+1, is generated by a group G e 02*+i with matrices gj e G 
have a form: 
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gi = [R l (i),...,R k (r>,±llV l<i<M. 

This implies that, in a code C(M, 2k + 1), the order M must be even, and the code is a union of two 
C(y,2fc) contained in parallel hyperplanes. Thus, an optimum commutative group code of order 
M in R 2k+l can be determined starting from the known optimal code in the previous dimension, 
C(f,2£), with initial vector 

xq = (S\,0,..., 6k, 0). The search for the best initial vector yg = (cos 9xq, sin 6) is then reduced 
to a single-parameter optimizing problem. 

The next section is devoted to the development of a method for searching for an optimum 
commutative group code free of reflection blocks in even dimensions. Besides providing addi- 
tional reduction in the number of cases to be tested, we show how to select and efficiently store 
a set of cases which allows to finding of an optimal code by solving the correspondent initial 
vector problems. 

4.1. Describing non-isometric commutative group codes 

Our approach starts with the connection between commutative group codes and lattices ll2~D . 
Specifically, let C be a commutative group code in R 2k , generated by a group G e B n j «, with 
matrices are free of 2 x 2 reflection blocks. We define the associated lattice Ac by 

Ac := [(b u . ..,b k )eZ k : [R(b,), . . .,R(b k )] e g) , 

where R(b) denotes the rotation in M. 2 by an angle of 2nb/M and [R(b\), . . . ,R(b k )] denotes a 
pseudo-diagonal matrix, according to ([TJ. 

We point out that Ac contains Ml} := {M(z\,zi, ■ ■ ■ ,Zk),Zi e Z) as a sub-lattice. Inside the 
hyperbox [0, M) k there are exactly M points of A G , which 
correspond to representatives of the elements of G, i.e., 

[0, M) k D {(b n ,b iz ,..., b ik ) mod M, i = 1 , 2, . . . , M) . 

The lattice Ac can then be viewed as the translation of these representatives through the 
lattice MZ k . 

If xq = (6\ , 0, . . . , 6k, 0) is an initial vector for the code C, we can also define a lattice A G (xo) 

by 
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A G (x ) := • 



b : b e A c 



Under these conditions, the code C is the image i^ A „(A G (xo)) c S 2k l , where 



t//,„(v) = \8\ cos(^-j,c5i sinf^- 



..ft cob VfWm | 



(2) 



is the standard parametrization of the torus with radii 5, ETI . 

We say that two lattices Aq and Ah are equivalent, and denote by Ac ~ Ah iff ift(Ac(xo)) 
and if/(A H (yo)) are isometric codes, for some xq, y e S 2k ~ l . 

As a consequence of the relation ~, we proceed to use isometry to discard commutative group 
codes to be checked in the searching for an optimum one. This will be done in terms of basis 
reduction of the associated lattices, based on results derived in Theorems |4.1|and|4.2| 



Theorem 4. 1 is closely related to a classical Hermite result. In particular, we have shown that 
the columns of the resulting matrix T can be ordered by the gcd (greatest common divisor) of 
their elements. In Theorem |4.2| we show that it is sufficient to consider generator matrices of 
lattices in a specific triangular form. 

Let Mi(%) be the set of k x k matrices with integer elements. GLi(Z) c M^Z) is the group 
of those matrices which are invertible in M,t(Z), the so called unimodular matrices. 

Theorem 4.1 (Special Hermite Normal Form). Let B be a k x k matrix with elements in Z. 
Then there is an upper triangular matrix T — U BV, with U € GLt(Z) and V a permutation 
matrix. Moreover, T satisfies the following conditions: 

1. < T(i, i) ^ T(i + 1, i + 1), V 1 < i < k - 1; 

2. sC T([l : i - 1], < T(i, /), V2<i<l; 

3. T(i, i) < gcd (T([i : j], J)) , V 1 < i < j < k; 

where T([p : q],r) are the elements in the rows p to q of the r-th column ofT. 
Proof: 

The proof is made by induction on k. For k — 1 it is trivial. Suppose the statement is valid 
for n < k. 
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Let Vi be a matrix which permutes the columns of B, such that the gcd of the column elements 
of the matrix B V\ are in increasing order. 

Let d\ = gcd((B V\\i) be the gcd of the elements in the first column of B V\, and U\ be a 
unimodular matrix, such that 



(3) 





d\ 










U\BVi = 




Bk,k-\ 




a k d\ 





i.e., the product of its first row by the first column of B V\ is equal to d\ . 
Let 



1 


••• 


-a 2 








-a k 





(4) 



be the matrix which provides the Gaussian elimination in the first column of BV 1 . We thus obtain 



f/if/i BVi = 



d\ 



















(5) 



Let B\ be the (k- 1) X {k - 1) submatrix of U\BV\, obtained by removing the first row and 

first column. By the induction hypothesis there exists a unimodular matrix U and a permutation 

matrix V such that T = U B\ V. 
Then, 
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1 




1 




U l BV l 




u 




V 



1 

u 



J, 







ch 


ch 




UBi 

eV 
UBi V 

eV 
f 



r. 



(6) 



If T(\,j) < or T(l,j) > T(j,j) for some j > 1, we can apply the elementary operation 



proof. 



U 



1 j-j 



(j, or equivalently left-multiply T by a unimodular matrix Uj, to conclude the 



In contrast to the standard Hermite normal form [35], here the unimodular matrix U is op- 
erating on the left side of B. In other words, if the rows of B contain the generator vectors of a 
^-dimensional lattice, then matrix U represents a change of basis in this lattice. Moreover, the 
permutation matrix V, which does not appear in the standard Hermite normal form, allows us to 
sort the columns of T by their greatest common divisor, which will be useful in order to discard 
isometric codes. We remark that the matrix V, operating on the right side of B, represents an 
isometry by coordinate permutation. Thus the lattices generated by T and B can be different, but 
they are isometric. 

Theorem 4.2. Every commutative group code C c S 2k ~ l , generated by a group G e 02k free of 
2x2 reflection blocks is isometric to a code obtained as image by tff of a lattice Ac(xo). Moreover 
the associated lattice Kg has a generator matrix T satisfying the following conditions: 

1. T is upper triangular according to Theorem \4.1\ 

2. de«T) = M k ~ l ; 

3. There is a matrix W, with integer elements satisfying WT — M 1^, where 1^ is the k x k 
identity matrix; 

M 

4. The elements of the diagonal of T satisfy T(i, i) — — where a,- is a divisor of M and 

ii, 

(fl;) ! • (fl, + i • ■ • Clk) < M, Vz* = 1, . . . , k. 
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Proof: 



1 - Let B be a generator matrix of the lattice A^. By Theorem 4.1 there exists an upper 



triangular matrix T such that T — U BV. Since the matrix U is unimodular, it defines a change 
of basis in the lattice generated by B, while V is an isometry by coordinate permutation. Both 
operations are isometries in lattices, thus, matrices B and T — U BV define lattices which are 
equivalent and which, therefore, generate isometric commutative group codes. 



2 - The lattice Ac contains the sublattice Ml k and the cardinality of the quotient 



A G 

Ml k 



must 



be equal to M, the number of points in the code. Therefore, since det(M4) = M k we conclude 
thatdet(r) =M k ~ l . 

3 - The system xT — M <?,- must have a solution in Z* for all 1 ^ i ^ k, where e, is the z'-th 
column of /j. Let W be the matrix with rows containing these solutions; then WT - MI^. ^ 

4 - The number M must be a multiple of the elements in the diagonal of T (from item 2) 
moreover, from Theorem |4T1 we know that 



M M 

T(i, i) ^ T(i + 1) and then — ^ which implies that a, + , ^ a,. 

a i cii+i 



From 



MM M , , 
de«T) = ... — = M , 



we get 



(«! 02 ■ • ■ at) — M => (a,)' • (a 



■ a k ) < M. 



Not all upper triangular integer matrices T satisfy the conditions of Theorem 4.2 For exam- 
ple, for M = 12 and k = 3, the matrix 



2 3 
6 6 
12 



Note that condition 3 is equivalent to saying that (Ml}) is a sublattice of A<j. 

12 



satisfies the hypothesis of Theorem 4. 1 and det(r) = 12 2 but, in order to obtain WT = 12 73, we 
must have: 



W 



6 -3 3/2 
2-1 
1 

However, in this case, W has non-integer elements. 

In order to characterize a commutative group code as an image of a quotient of lattices, it 
is also important to determine a set of generators of the correspondent group and its class of 



isomorphism. In the Theorem 4.5 we deal with this problem. 



Theorem 4.3 (|35|, p 76). Let A be a non-singular k x k matrix with coefficients in Z. There is 
then a unique diagonal matrix D — (djj), with di+ij+\\dij, such that D — V AU with U and V in 
GL k (T). 

This matrix is called the Smith normal form (SNF) of A. 

Theorem 4.4 (|35|, p 76). Let L be a 'Z-submodule of a free module L' and of the same rank. 
Then there are positive integers d\, . . . ,c4 satisfying the following conditions: 

1. For every i such that 1 ^ i < k we have di + \ \di. 

2. As ^-modules, we have the isomorphism 



and in particular [L' : L] = d\- • • du and d\ is the exponent of 4-. 
3. There is a Z-basis fvi, . . . , Vjt} ofL' such that {d\y\,..., dk\k) is a 2,-basis ofL. 

Theorem 4.5. For a commutative group code C, let T be a generator matrix of the lattice Ac, 
according Theorem 4.2 and W = MT^ 1 . The set of generators of the correspondent group and 



its class of isomorphism are then obtained from the SNF ofW. 



Proof: 

Let D = V W U be the SNF of W. We know that 
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WT = MI k => V~ l DU l = MI k => D U~ l T = VMI k . 

Since the matrices U~ l e V are unimodular, their product on the left of the generator matrices 
T and M I k define a change of basis in the lattice generated by T and its sublattice Ml}. The 



classification and generators of the group are derived from Theorem 4.4 In this case, G is 
isomorphic to a group Z^, ® . . . ® Z</ t and the rows of I/ -1 T give the elements bjj which form a 
set of generators, according to ([TJ. 



As a consequence of these results, we derive a two-step algorithm which searches for an 
optimum commutative group code C of order M in an even dimension. The first step consists 
of storing all matrices T according to theorem 



4.2 



and the use of Adam's relation to discard 



isometric groups. For each one of these matrices T we then establish a linear programming 
problem (Section 3) to determine the initial vector Xq which maximizes the minimum distance 
of the group code i!j Xo Aq{xq) Q. For the optimum case, theorem 4.5 is applied to obtain the 
generators and the class of isomorphism of the commutative group. The algorithm is summarized 
as a pseudo code in Algorithm 1. 

Let us illustrate this method in detail for M = 128 and n — 4. 

Let div = {1,2,4,8, 16,32,64, 128} be the set of divisors of 128. From Theorem |42j we 
know that the matrix T, related to a code C(128,4), has the form 



d\ w 
d 2 



■UJ M 

, with dj = — , fl; e div. 

at 



Moreover (c/2) 2 < 128, i.e., d 2 e {1,2,4,8}. We can then store all the possible diagonal of T as 
columns of a matrix A: 



A = 



12 4 8 
128 64 32 16 

For each column of A, the set of values w in T can then be determined, as established in item 3. 



of Theorem 4. 1 by considering 



a h sC gcd(w,a 2 ,), 
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Algorithm 1: Optimum commutative group code 



input : The number of points M and the dimension n = 2k; 

output: An optimum commutative group code C(M, n), its set of generators, optimum 
initial vector, its isomorphism class and minimum distance. 



begin 

dist 



0; 



div <— «2, . . . , a w ], the set of divisors of M; 

A <— \diagi\, diagi\, c//ag ; |], a matrix with columns contains all the possible 
diagonals for T, according Theorem |4. 2 1 i.e., 

diagi = ^, . . . , ■ffi, where a uk e div, a iJc > a iMl and \\ k 

foreach diagi € A do 

Step 1: Construct all matrices T, according Theorem 

to discard isometric groups; 



4.2 



?=1 a Uq = M; 



and use Adam's relation 



foreach matrix T£ constructed in step 1 do 

Step 2: Solve the initial vector problem and get the minimum distance dist^ 

and the initial vector xq', 

if distil > dist then 
dist <— distjf, 



T < 



Apply Theorem 4.5 in T and get the generator of the group G e O n and the 

correspondent isomorphism class; 

Output G, xq, dist and the isomorphism class. 
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In this example, we have 



r, = 



Ta = 



1 







128 


2 







64 


4 







32 


8 







16 



where w\% e {0, 1, • • • ,64} ; 
where w n e {0, 2, 4, 6, • • • , 32} ; 
where w 3f e {0,4,8, 12, 16); 
where e {0, 8} . 



This amounts to 89 cases to be tested. However some of these lattices are equivalent. For 
example, the lattice generated by a matrix T\ which has the first row equal to (1, w\^) is equivalent 
to a lattice generated by a matrix T\ which has the first row equal to (1, Wji), here W\l represents 
the inverse of w\% in 7L\ 2 %. If < w\%, we can therefore discard the correspondent matrix in set 
T\ . This situation occurs for 



w X( = {17,27,33,35,39,41,43,45,49,51,53,55,57,59,61}. 



Similarly, the lattice generated by a matrix T 2 , which has the first row equal to (2, 2b) is 
equivalent to a lattice generated by a matrix T 2 which has the first row equal to (2, 2b~ l ). Thus, 
in this set we can discard the cases where = {18, 26). Therefore, in order to find an optimum 
code C(128, 4) it is sufficient to check 72 codes. 

In the implementation of Algorithm 1, these equivalent cases can be discarded during Step 
1 and the solution of the initial vector problem, consequently, implemented just for the relevant 
cases. Only the matrix which determines the largest minimum distance must be saved. 

1 11 

128 

The correspondent group G e <9 n is then obtained using the SNF of 
W = M(T U2 y\ 

In this case, the best C(128,4) is a cyclic group code with the following generator matrix: 



In this example, the optimum code is associated to the matrix T\xi 
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r \*2tt\ . l*2n\ 1 

cos sin 

\ 128 ) \ 128 ) 

(l*2n\ (l*2n\ 

- sin cos 

r - V 128 / V 128 / 

0(1,11,128) - /ll»2w\ /ll*2w\ • 

cos sin 

V 128 / \ 128 j 

. Ill*2n\ {ll*2n\ 

- sin cos 

L \ 128 ) \ 128 ) . 

The minimum distance in this code is d — 0.406179 for the best initial vector 
x Q = (0.65098,0,0.759095,0)'. 

In dimension 4, the number of commutative group codes tested by Algorithm 1 is not much 
larger than the number of cyclic group codes tested (Table [TJ. For M equals to 32 (respectively, 
64,128,256,512,1024), Algorithm 1 checks 21 (respectively, 38,72,141,273,542) commuta- 
tive group codes in order to find an optimum one. 

Using this method we have found optimum codes for various values of M in different dimen- 
sions and we have present some of them in M 4 and R 6 in Tables|2]and[3] In both cases, it can be 
seen that, when the number of points M increases, the gap between the minimal distance of the 
codes and the upper bound ll2Tl decreases. This fact is also illustrated in Figure [T] 



Table 2: Some optimum commutative group codes of order M in R 4 . 



M 




Si 


Si 


Group 


Gen. (bij) 


Bound 


10 


1.224 


0.707 


0.707 


Zio 


(13) 


1.474 


20 


0.959 


0.678 


0.734 


Z20 


(3 4) 


1.054 


30 


0.831 


0.707 


0.707 


Z30 


(3,5) 


0.864 


40 


0.714 


0.607 


0.794 


Z40 


(4 5) 


0.750 


50 


0.628 


0.707 


0.706 


Z50 


(7 2) 


0.672 


100 


0.468 


0.757 


0.653 


z 5 e I20 


(0 20), (5 10) 


0.476 


200 


0.330 


0.750 


0.660 


Z200 


(93 1) 


0.337 


300 


0.273 


0.656 


0.754 


z 5 © i 60 


(60 120), (10 15) 


0.275 


400 


0.237 


0.686 


0.727 


Z400 


(189 1) 


0.238 


500 


0.211 


0.674 


0.738 


Z500 


(13 20) 


0.213 


600 


0.193 


0.676 


0.736 


^600 


(191 198) 


0.194 


700 


0.180 


0.718 


0.695 


^700 


(14 25) 


0.180 


800 


0.168 


0.670 


0.742 


^800 


(16 25) 


0.168 


900 


0.158 


0.704 


0.709 


^900 


(197 2) 


0.159 


1000 


0.149 


0.716 


0.697 


^1000 


(33 4) 


0.150 
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Table 3: Some optimum commutative group codes of order M in R 6 . 



M 


d m ,„ 




Si 


<S 3 


Group 


Gen (bjj) 


Bound 


10 


1.414 


0.632 


0.632 


0.447 


Zio 


(3,1,5) 


1.820 


20 


1.240 


0.554 


0.620 


0.554 


Z 20 


(2,5,6) 


1.465 


30 


1.133 


0.534 


0.654 


0.534 


Z 3 o 


(3,5, 9) 


1.287 


40 


1.044 


0.603 


0.522 


0.603 


z 2 e z 20 


(20,0,20), (32,10,4) 


1.173 


50 


0.976 


0.604 


0.506 


0.615 


z 50 


(7,6, 34) 


1.091 


100 


0.804 


0.515 


0.684 


0.515 


z, © Zio 


(50, 10, 0), (30,0, 10) 


0.870 


200 


0.673 


0.555 


0.619 


0.555 


z 2 oo 


(28, 25, 4) 


0.692 


300 


0.585 


0.585 


0.498 


0.639 


Zs ffi z m 


(0, 0, 60), (25, 30, 30) 


0.605 


400 


0.540 


0.562 


0.605 


0.562 


z 20 © z 20 


(300, 40, 0), (60, 0, 20) 


0.550 


500 


0.504 


0.577 


0.577 


0.577 


Z 5 © Zm, <8Z10 


(100. 0, 0), (50, 50, 0), (50, 0, 50) 


0.511 


600 


0.472 


0.549 


0.630 


0.549 


z 2 © Z300 


(300, 0, 300), (384, 50, 12) 


0.481 


700 


0.445 


0.531 


0.612 


0.585 


Z700 


(457, 664, 298) 


0.457 


800 


0.427 


0.617 


0.486 


0.617 


z 2 o © z 40 


(80,0,40),(20,80,60) 


0.437 


900 


0.413 


0.592 


0.591 


0.547 


Z 3 ffl z 300 


(0,300,0),(759,36,3) 


0.420 


1000 


0.397 


0.560 


0.632 


0.535 


Z1000 


(319,694,45) 


0.406 



Bound for d mln in R* 




50 100 150 200 250 300 

Figure 1: Comparison between the distance of optimal codes found using Algorithm 1 and upper bound |21|: the gap 
decreases when M grows. 

Although some other group codes, as permutations codes, can outperform commutative 
group codes for some parameters fl9l . they are very special for some applications as transmition 
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over symmetric channels IfTTl . Besides, they may provide homogeneous spherical codes for any 
number of codewords and can be used for designing high density spherical codes on fiat torus 
layers (29). 

5. Conclusions 

A two-step method for finding an optimum n-dimensional commutative group code of order 
M is presented. The approach explores the structure of lattices associated with these codes in 
even dimensions and allows a significant reduction in the number of non-isometric cases to be 
analyzed. For each of these cases, a linear programming problem is solved to find the initial 
vector which maximizes the minimum distance in the code. The method introduced here can 
also be used to design more general spherical codes, such as the so called quasi-commutative 
group codes, which are constructed on layers of flat tori [29 1. 
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